Spin wave excitation in magnetic insulators by spin-transfer torque 
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We study the excitation of spin waves in magnetic insulators by the current-induced spin-transfer 
torque. We predict preferential excitation of surface spin waves induced by an easy-axis sur- 
face anisotropy with critical current inversely proportional to the penetration depth and surface 
anisotropy. The surface modes strongly reduce the critical current and enhance the excitation 
power of the current-induced magnetization dynamics. 



Spintronics is all about manipulation and transport of 
the spin, the intrinsic angular momentum of the electron 
These two tasks are incompatible, since manipula- 
tion requires strong coupling of the spin with the outside 
world, which perturbs transport over long distances. In 
normal metals spin can be injected and read out eas- 
ily, but the spin information is lost over short distances 
[2J. In spin-based interconnects, transporting spins over 
longer distances is highly desirable Q. 

The long-range transport of spin information can be 
achieved by encoding the information into spin waves 
that are known to propagate coherently over centime- 
ters [4j. It has been demonstrated in Refs. MtA for the 
magnetic insulator Yttrium-Iron-Garnet (YIG) that spin 
waves can be actuated electrically by the spin-transfer 
torque [H 0| and detected by spin pumping (1(1 at a dis- 
tant contact. In the experiment by Kajiwara [5j, Pt was 
used as spin current injector and detector, making use 
of the (inverse) spin Hall effect 11] . In a d = 1.3 /im- 



thick YIG film spin waves were excited by a threshold 
charge current of J c ~10 9 A/m 2 . This value is much less 
than expected for the bulk excitation that in a linear 
approximation corresponds to the macrospin mode and 
is estimated as J c = (l/6 H )eauM s d/jh~lO n ~ u A/m 2 , 
where e and 7 are the electron charge and gyromagnetic 
ratio, respectively, and we used the parameter values in 
Table U for the ferromagnetic resonance frequency u>, the 
spin Hall angle of Pt 9h, magnetic damping a, and sat- 
uration magnetization M s . 

In this Letter, we address this large mismatch be- 
tween observed and expected critical currents by study- 
ing the threshold current and excitation power of current- 
induced spin wave excitations. We present a possible an- 
swer to the conundrum by proving that the threshold 
current is strongly decreased in the presence of an easy- 
axis surface anisotropy (EASA). Simultaneously, EASA 
increases the power of the spin wave excitation by at least 
two orders of magnitude. 

We study a structure as depicted in Fig. [TJ where a 



non-magnetic (N) metallic thin film of thickness t is in 
contact with a ferromagnetic insulator (FI), whose equi- 
librium magnetization is along the ^-direction. The spin 
current injected into the ferromagnetic insulator is polar- 
ized transverse to the magnetization J s = J s mxzxm. 
The bulk magnetization is described by the Landau- 
Lifshitz-Gilbert (LLG) equation: 

rh = -7DIX [H + (Aex/i) V 2 m + h] + amxm, (1) 

where Ho includes the external and internal magnetic 
field, A cx is the exchange constant, and h is the dipo- 
lar field that satisfies Maxwell's equations. In the qua- 
sistatic approximation, i.e. disregarding retardation in 
the electromagnetic waves, Vxh = and V-b = V-(h + 
fXoM s m) — 0. All quantities are position and time depen- 
dent. In the absence of pinning, the total torque vanishes 
at the interface (l2j : 
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where n is the outward normal as seen from the ferro- 
magnet. The first term in Eq. ([2]) is the surface exchange 
torque, the second term the torque due to a perpendic- 
ular uniaxial surface anisotropy H a = ^p(m-n)n and 
K s = JdxKi across the surface, and the last term is the 
current- induced spin-transfer torque [13j . We parame- 
terize the surface anisotropy and spin current as wave 
numbers k s — 2"fK s / A ex M s and kj — 7J S / A CX M S . The 
dipolar fields h y , z and b x are continuous across the inter- 
face. Eqs. ((lj [2j) in combination with Maxwell's equations 
describe the low energy magnetization dynamics and can 
be transformed into a 6th-order differential equation for 
the scalar potential ip with h 



VV> PJ,ll5|. 



The method described above extends a previous study 
by Hillebrands 15] by including the current-induced spin- 
transfer torque. We predict the critical conditions under 
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FIG. 1. (Color online) An electrically insulating magnetic film 
of thickness d with magnetization m (||z at equilibrium) in 
contact with a normal metal. A spin current J s ||z is generated 
in the normal metal and absorbed by the ferromagnet. 

which magnetization dynamics becomes amplified by the 
current-induced driving torque. 

We start with the limiting case of d — > oo (semi-infinite 
ferromagnet). After linearization and Fourier transfor- 
mation in both time and space domains, Eq. ([1} reduces 
to a 4th-order differential equation in ip. Focusing for 
simplicity first on the case of vanishing in-plane wave- 
vector q = (q y , q z ) = 0, the scalar potential can be writ- 
ten as: V( r ) = J }2j=i a j e%qjX e%ult with 



<?jH = -i 



UJ 



w 2 + \oj 2 



M 



(3) 



and | gi | ^192 1 when u^loq. Imposing the boundary con- 
dition in Eq. ([2]) , up to the first order in kj : 



= 2q 1 q 2 (qi + q^)+ik s 



{qi + q 2 ) 



■4k jU . (4) 



The solutions of Eq. (U) are the complex eigen- frequencies 
to, whose real part represents the energy and imaginary 
part the inverse lifetime. To Oth-ordcr in dissipation, i.e. 
with vanishing bulk damping (a = 0) and spin current 
injection (kj — 0), and using |gi|^|(?2|, Eq. Q simplifies 
to k s — iq 2 / [l+u>M / (Aql)], which has no non-trivial solu- 
tion for k s < 0. The single real solution for k s > obeys 
uj < \J loq(ujo + wjj) such that both q\. 2 are negative 
imaginary: q X — ~ i \/(2ujq + w M )/A cx , q 2 ~-ik s uj Q /(2u! + 
+ O(k^), i.e. a surface spin wave induced by the 
easy-axis surface anisotropy. With the criteria Im u) < 
and to leading order in < k s <^qi, Eq. ([4]) leads to the 
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YIG 


Unit 
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1.76xlO n 


1/(T s) 


M s 


a 1.56xl0 5 


A/m 


u« = J Ho M s 


34.5 


GHz 


Aex 


4.74xl0" 6 


m 2 /s 


a 


a 6.7xl0 -5 




wo = jHo 


0.5wm 


GHz 


Ks 


b 5xl(T 5 


J/m 2 



TABLE I. Parameters for YIG. a RefJE b K s ranges 0.01~0.1 
erg/cm 2 or 10~ 5 ~10 -4 J/m 2 , Ref. 0Jand[l7|. 



critical current: 

a (w +w M /2) 2 
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When there is no surface anisotropy (k s — > 0), the criti- 
cal current diverges because the macrospin mode cannot 
be excited in a semi-infinite film. Using the parameters 
given in TableUin Eq. ([5]), we estimate the critical current 
for exciting the EASA induced surface wave (at q = 0) 
to be kj = — 0.08A; C , where k c = a(ujQ + uim /2)d/A e * is 
the critical current for bulk excitation in a YIG thin film 
of thickness d = 0.61 /j,m (used below). 

EASA pulls down a surface spin wave for the following 
reason: when k 3 ■ — J s = 0, the boundary condition in 
Eq. @ requires cancellation between the exchange and 
surface anisotropy torques: d x m x — k s m x — d x m y — 
0. The exchange torque depends on the magnetization 
derivative in the normal direction, and can only take one 



sign in the whole film, and m x 



as x 



there- 



fore (l/m x )d x m x > 0. Torque cancellation (for a non- 
trivial solution) is therefore possible only for k s > 0. The 
surface spin wave induced by EASA (k s > 0) for the in- 
plane magnetized film (m z ~l) discussed in this Letter 
is analogous to the surface spin waves for the perpen- 
dicular magnetized film (m x ~\) induced by easy-plane 
surface anisotropy (k s < 0) studied before in YIG films 
[l2l Il8l - [2"H . For perpendicular magnetization, a different 
boundary condition: d x m y , z + k s m y ^ z = results in a 
surface wave for k s < 0. 

We now include all ingredients: finite thickness (d = 
0.61 //m), surface anisotropy, intrinsic magnetic damp- 
ing, spin current injection, exchange coupling, and dipo- 
lar fields. We calculate numerically the complex eigen- 
frequencies w(q, kj) as a function of the in-plane wave- 
vector q and the applied spin current at the surface kj. 
Im ui, the effective dissipation, can be either positive 
(damping) or negative (amplification) when driven by the 
spin-transfer torque. 

First, we disregard the surface anisotropy: K s — k s — 
0. With 9 the angle between q and m, the results for 
9 = 90° are shown in Fig. [2l In the top left panels Re u, 
the magnetostatic surface wave (MSW) is seen to cross 
the flat bulk bands [3] ■ When no spin current is applied 
(kj = 0), the dissipative part Im w^a(wo + Wm/2) > 0, 
as shown in the top middle panels. At a spin current that 
is 20% of that required for bulk excitation: kj — 0.2fc c , 
the dissipative part Im u (top right panel) decreases while 
Re u remains unchanged because the spin-transfer torque 
as magnetic (anti-)damping mainly affects Im w. Nega- 
tive effective dissipation implies spin wave amplification. 
This happens for the 5th (green) band at gde[2,6.5], 
which corresponds to a (chiral) MSW (mixed with bulk 
modes) formed near the interface (shaded/yellow panel). 
On the other hand, for 9 — —90° (not shown), the mag- 
netostatic surface wave at the opposite surface to vacuum 
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FIG. 2. (Color online) Spin wave band structure and mag- 
netization profiles in YIG for d — 0.61 /J,m without surface 
anisotropy: k s — at 6 — Z(m,q) = 90°. Top (from left to 
right): Re(w/wjw) vs. qd, Im (w/Vm ) at fcj = 0, Ihi(w/wm) 
at fcj = — 0.2fc c . Bottom: of the same 6 modes for 

qd — 0.09 and 3.74 indicated by the dashed vertical lines in the 
top panels. The colors label different bands. The amplitude 
of the green mode mode (shaded/yellow panel) is amplified. 



(x = — d) is only weakly affected by the spin current in- 
jection at x = 0. 

We now turn on EASA: k s = 25.0/^m (or K s = 
5xlCT 5 J/m 2 ) at the top surface (x = 0). Fig. [3] shows 
the results for = 90°. The changes of Re lu and Im lu 
at kj — are modest (Fig. [5]), but an additional band 
(black) appears, viz. the surface spin wave band in- 
duced by EASA. The spin-transfer torque strongly af- 
fects this mode because of its strong surface localization 
22|. As seen in the top right panel, almost the whole 
band is strongly amplified by a spin current injection of 
kj = 0.2k c . Inspecting the spin wave profiles at two dif- 
ferent q values, we observe a surface spin wave near x = 
for the black band at small q (shaded/yellow panel in the 
middle row in Fig. [3]). At larger q, the 1st (black) band 
loses its surface wave features to the 5th (red) band (see 
top right panel in Fig. [3]). The red band mode starts 
out as a magnetostatic surface spin wave, but the EASA 
enhances its surface localization by hybridization with 
the black mode to become strongly amplified by the spin 
current at higher q. Also in the lower panel of Fig. [3] 
we observe that the red band has acquired the surface 
character. 

We introduce an approximate power spectrum (Fig. |4]) 
that summarizes all information about the mode- 
dependent current-induced amplification: 



P(w) = 



Im w„<0 



|Im w n (q)|<5[w-Re w„(q)]dq (6) 
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FIG. 3. (Color online) Same as Fig. [2]but with k s — 25//im. 



with n the band index is the density of states at fre- 
quency lu weighted by its amplification. Without sur- 
face anisotropy, only a few modes are excited even at a 
relatively large current (kj — 0.2k c ). However, when 
k s = 25//Ltm, the excitation is strongly enhanced by 
more than two orders of magnitude due to the easily 
excitable surface spin wave modes. Furthermore, we ob- 
serve broadband excitation over a much larger range of 
frequencies. This power spectrum is rather smooth, while 
the experiments by Kajiwara et al. [5j show a large num- 
ber of closely spaced peaks. The latter fine structure is 
caused by size quantization of spin waves due to the finite 
lateral extension of the sample that has not been taken 
into account in our theory since it complicates the calcu- 
lations without introducing new physics. The envelope 
of the experimental power spectrum compares favorably 
with the present model calculations. 
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FIG. 4. (Color online) Top: Power spectrum (resolution 
8u/ujm = 0.01) at various current levels (kj = 0.2k c from 
the top decreasing by Akj — 0.01fc c ) without (left: k s = 0) 
and with (right: k s — 25.0/^tm) surface anisotropy. Inset: the 
integrated power versus kj. 
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The insets in Fig. @] show the integrated power and 
allow the following conclusions: 1) the excitation power 
is enhanced by at least two orders of magnitude by the 
EASA; 2) the critical current for magnetization dynamics 
is fcj^0.08fc c for k s = 25/^m, which agrees very well with 
the estimates from Eq. (|S|). This critical current is about 
one order of magnitude smaller than that for the bulk 
excitation (k c ), and about half of that for MSW with- 
out surface anisotropy (kj = 0.16/c c ). For k s = 25//im, 
it corresponds to J c = 3xl0 10 A/m 2 for 6 H = 0.01 Q 
and 3.8xl0 9 A/m 2 for 9 H = 0.08 [Mill. Th ese values 
are calculated for a film thickness of d = 0.61 /im, but 
should not change much for d = 1.3 /jm corresponding 
to the experiment Q , because the excited spin waves are 
localized at the interface. Compared to the original esti- 
mate J c ^10 11 ^ jl2 A/m 2 , the critical current for a surface 
spin wave excitation is much closer to the experimental 
value of J C '~10 9 A/m 2 [f| (although these experiments re- 
port a very inefficient spin wave absorption in contrast 
to the present model assumption). 

According to Eq. ([5]), critical current (excitation 
power) would be further reduced (increased) by a larger 



EASA. Ref. Lly reports an enhancement of the YIG sur- 



face anisotropics for capped as compared to free surfaces. 
A Pt cover on a YIG surface [f| may enhance the sur- 
face anisotropy as well. As seen from Fig. [3j the surface 
mode (black band) has group velocity duj/dq compara- 
ble to that of the MSW. The excited surface spin wave 
therefore propagate and can be used to transmit spin in- 
formation over long distance at a much lower energy cost 
than the bulk spin waves. 

In conclusion, we predict that an easy-axis surface 
anisotropy gives rises to a surface spin wave mode, which 
reduces the threshold current required to excite the spin 
waves and dramatically increases the excitation power. 
Multiple spin wave modes can be excited simultaneously 
at different frequencies and wave- vectors, thereby ex- 
plaining recent experiments. Surface spin wave excita- 
tions could be useful in low-power future spintronics- 
magnonics hybrid circuits. 
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